Spin waves in a one-dimensional spinor Bose gas 
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We study a one-dimensional (iso)spin 1/2 Bose gas with repulsive (5- function interaction by the 
Bethe Ansatz method and discuss the excitations above the polarized ground state. In addition 
to phonons the system features spin waves with a quadratic dispersion. We compute analytically 
and numerically the effective mass of the spin wave and show that the spin transport is greatly 
suppressed in the strong coupling regime, where the isospin-density (or "spin-charge") separation 
is maximal. Using a hydrodynamic approach, we study spin excitations in a harmonically trapped 
system and discuss prospects for future studies of two-component ultracold atomic gases. 

PACS numbers: 



Recent experiments have shown the possibility of 
studying ultra-cold atomic gases confined in very elon- 
gated traps In such geometries, the gas be- 
haves kinematically as if it were truly one-dimensional 
(ID). Many theoretical studies H IS II H II 113 have 
predicted and discussed interesting effects in ID Bose 
gases, such as the occurrence of fermionization in the 
strong coupling Tonks-Girardeau (TG) regime, where el- 
ementary excitations are expected to be similar to those 
of a non-interacting ID Fermi gas (5|- Manifestations of 
strong interactions have been found in the experiments 
0, and recently the TG regime has been achieved for 
bosons in an optical lattice Q and in the gas phase 0. 

Present facilities allow one to create spinor Bose gases 
which has been demonstrated in experimental studies of 
two-component Bose-Einstein condensates 0|. These 
systems are produced by simultaneously trapping atoms 
in two internal states, which can be referred to as 
(iso)spin 1/2 states. Relative spatial oscillations of the 
two components can be viewed as spin waves |l2j (see 
|l3| for review) . A variety of interesting spin-related ef- 
fects such as phase separation 1^, exotic ground states 
|l5|. and counter intuitive spin dynamics 12] due to the 
exchange mean field, have been studied both theoreti- 
cally and experimentally. However, most of these studies 
are restricted to the weakly interacting Gross-Pitaevskii 
(GP) regime. There is a fundamental question to what 
extent these effects survive in the strongly correlated 
regime characteristic of one spatial dimension. The pur- 
pose of the present Letter is to study spin excitations of 
an interacting ID spinor Bose gas. This is done using an 
exact solution by the Bethe Ansatz. 

We start with a spinor (two component) gas of N 
bosons with mass m at zero temperature, interacting 
with each other via a repulsive short-range potential 
in a narrow three-dimensional waveguide. In general, 
the interaction depends on the internal (spin) states of 
the colliding particles. Here we consider the case of a 



spin-independent interaction characterized by a single 3D 
scattering length a > 0. This is a reasonable approxima- 
tion for the commonly used internal levels of *^Rb (see 
e.g. The waveguide has length i and we assume pe- 

riodic boundary conditions for simplicity. The transverse 
confinement is due to a harmonic trapping potential of 
frequency ujq. When the chemical potential of the gas is 
much smaller than hujQ , the transverse motion is frozen to 
zero point oscillations with amplitude Iq = yJh/mujQ. In 
such a quasi-lD geometry, the interaction between atoms 
is characterized by an effective one-dimensional delta- 
potential g5(x). For a <C ^O; the coupling constant g is 
related to the 3D scattering length as g = 2h^a/mlQ > 
. The behavior of the system depends crucially on the 
dimensionless parameter 7 = mg/h^n, where n = N/L 
is the ID density. For 7 <C 1 one has the weak coupling 
GP regime, whereas for 7^1 the gas enters the strongly 
interacting TG regime. 

Under the above conditions, the system is governed by 
the following spin-independent ID total Hamiltonian: 
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This Hamiltonian was introduced by Lieb and Liniger 
for describing spinless bosons, and their solution by 
the Bethe Ansatz (BA) has been generalized to bosons 
or fermions in two internal states by M. Gaudin and 
C.N. Yang |lHlr^ . In the case of a two-component Bose 
gas (spin 1/2 bosons), due to the SU{2) symmetry of 
the Hamiltonian the eigenstates are classified according 
to their total (iso)spin S ranging from to N /2. In this 
case, which was recently considered by Li, Gu, Yang and 
Eckern J^], the ground state is fully polarized [S — N/2) 
and has 25* -|- 1-fold degeneracy, in agreement with a gen- 
eral theorem 0,113. At a fixed S = N/2, the system is 
described by the Lieb-Liniger (LL) model 0, for which 
elementary excitations have been studied in Ref. 0] for 
any value of the interaction constant. Spin excitations 



2 



above the ground state are independent of the ground- 
state spin projection Ms and represent transverse spin 
waves. For Ms = they correspond to relative oscilla- 
tions of the two gas components. 

We fir st g ive a brief summary of the BA diagonaliza- 
tion 0, 0i 01 of the Hamiltonian An eigenstate 
with total spin S = N/2 ~ K (0 < K < N/2) is char- 
acterized by two sets of quantum numbers: N density 
quantum numbers Ij with j ~ I, N and K spin quan- 
tum numbers with /i = 1, .., K. If N — K is odd (resp. 
even), Ij and are integers (resp. half- integers) . These 
quantum numbers define N quasi-momenta kj and K 
spin rapidities , which satisfy the following set of BA 
equations (we set h = 2m = 1): 
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The energy of the corresponding state is E' = fc| , and 
its momentum is given by: 
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As we are also interested in finite size effects, we do not 
take the thermodynamic limit at this point. 

The ground state corresponds to the quantum numbers 
{/O} = {-(AT _ i)/2, .., {N - l)/2} and if = 0, which 
shows that the BA equations reduce to those of LL 0. 
The wave function is given by the orbital wave function of 
the LL ground state multiplied by a fully polarized spin 
wave function. All ground state orbital properties (en- 
ergy, chemical potential, correlation functions, etc.) are 
therefore identical to those of the LL model. Elementary 
excitations in the density sector correspond to modifying 
the density quantum numbers Ij while leaving the total 
spin unchanged, i.e., -ftT = 0. At low energy, the density 
excitations are phonons propagating with the Bogoliubov 
sound velocity Vg = \J2gn in the GP limit and with the 
Fermi velocity Vs = 27m in the TG regime. 

We now focus on the spin sector. Elementary spin exci- 
tations correspond to reversing one spin [K = 1), and the 
total spin changes from N/2 to A^/2 — 1. Thus, we have 
a single spin rapidity A and the corresponding quantum 
number J. In general, this procedure creates a density 
excitation and a spin wave (isospinon) jisj . Here, we 
choose specific quantum numbers Ij , J in order to excite 
the isospinon alone |2ll |. Accordingly, the momentum p 
of the excitation is 



27r f N 



(5) 



which follows from the definition Q). 



In the long wavelength limit, where \p\ <C n, due to 
the SU{2) symmetry one expects ^ a quadratic disper- 
sion for the spin-wave excitations above the ferromag- 
netic ground state: 



E{p) -Eo^ //2m* 



(6) 



where E{p) is the energy of the system in the presence of 
a spin wave with momentum p, Eq is the ground state en- 
ergy and m* is an effective mass (or inverse spin stiffness) . 
This quadratic behavior is due to a vanishing inverse 
spin susceptibility, which is a consequence of the SU (2) 
symmetry 22j. A variational calculation in the spirit of 
Feynman's single mode approximation |2Clj . shows that 
£p ^ p^ /2m implying that m* > m. Below we show that 
strong interactions greatly enhance the effective mass. 

In the strong coupling limit it is possible to solve the 
BA equations Q and © perturbatively in I/7 |^. We 
solve these equations both for the ground state {I^} and 
the excited state {Ij] J}. We anticipate that in the limit 
of strong interactions, for small momenta (|p|/?T. «C 1) and 
a large number of particles {N 3> 1), the dimensionless 
spin rapidity is A = 2A/g 3> 1 and the dimensionless 
quasi-momenta arc \kj\/g <^ 1. This allows us to expand 
Eqs. (0) and Q to first order in I/7 and 1/iV. The 
ground state quasi-momenta are then given by: 



kU = 2nl" (1 - 2/7) . 



(7) 



Here we used the relation arctan(2(fc? — kf)/g) ~ 
2Nkj/g — 2{J2i k'^) / 9 = '^Nkj/g, which is a consequence 
of the vanishing ground state momentum. Similarly, the 
excited state quasi-momenta obey the equations: 
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where p is given by Eq. Q. Neglecting quasi-momenta 
ki in the argument of arctangent in the BA equation 
we obtain the excited state spin rapidity: 

27r J = 2N arctan(2A) ~ ttTV - N/\ (9) 

Equations ||SJ| and then give: 

A = N/pL, (10) 

which justifies that A ^ 1 for \p\/n <C 1. Combining this 
result with Eq. © shows that \kj\/g ^ 1, as anticipated. 
Let us now define the shift of the quasi-momenta Akj = 
kj — k'j. Taking the difference between equations (jSJ and 
0, we find: 
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77VA2 



2p_ 

jN 



2tt 
i7 



(11) 



1/2. We can now compute 



the energy of the spin wave, as defined in Eq. JBJ 
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Using Eq. (fTT|l for Afcj and Eq. lfTU|l for A gives £p = 
p2 (l/iV + 27r2/37). Note that the last two terms in the 
right hand side of Eq. give no contribution, as the 
ground state momentum is zero. According to the defi- 
nition ©, the inverse effective mass is therefore: 



m 
m* 



1 

N 



2n^ 
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where we restored the units. Remarkably, the effective 
mass reaches the total mass Nm for 7 00: the bosons 
are impenetrable and therefore a down spin boson can 
move on a ring only if all other bosons move as well. 

In the opposite limit of weak interactions it is possi- 
ble to conipute the effective mass from the Bogoliubov 
approach |23| . The validity of this procedure when con- 
sidering a ID Bose gas, i.e. in the absence of a true 
Bose-Einstein condensate, is justified in The Hamil- 
tonian of the system can be written as Hq + Hint, where 
Hq is the Hamiltonian of free Bogoliubov quasiparticles 
and free spin waves: 



Hn 



p 



(14) 



with ap, Pp being the Bogoliubov quasiparticle and the 
spin wave field operators, ep = ^ ep{ep + 2gn) the Bo- 
goliubov spectrum, and Cp — /2m the spectrum of 
free spin waves jl^ . The Hamiltonian Hint describes the 
interaction between Bogoliubov quasiparticles and spin 
waves and provides corrections to the dispersion relations 
tp and Cp. The most important part of Hint reads: 



Hint — 



PUf3k 
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where Uq and Vq are the u, v Bogoliubov coefficients sat 
isfying the relations Uq 



q — -\/ £q/ Cq and Uq — Uq — 

Neglected terms contribute only to higher 



orders in the coupling constant. To second order in per- 
turbation theory, in the thermodynamic limit the pres- 
ence of a spin wave changes the energy of the system by: 
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In order to calculate a correction to the effective mass of 
the spin wave, we expand Eq. (|16() in the limit of p — > 0. 
Terms which do not depend on p modify the ground state 
energy, linear terms vanish, and quadratic terms modify 
the spin wave spectrum as follows: 
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Ag^n 



dq 
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(17) 



where the main contribution to the integral comes from 
momenta q ~ y/mgn. Using the definition lO, we then 




FIG. 1: Inverse effective mass ra/m* as a function of the 
dimensionless coupling constant 7 (logarithmic scale). The 
stars (*) show numerical results for TV = 111 particles, the 
solid curve represents the behavior in the strong coupling limit 
(Eq. 113ll 'l. and the dashed curve the behavior for a weak 
coupling (Eq. CHJ)- 



obtain the inverse effective mass: 



m _ ^ 2^ 
m* TT 



37r 



which clearly shows non-analytical corrections to the bare 
mass due to correlations between particles. This result 
can also be obtained directly from the BA equations. 

For intermediate couplings, we obtained the effective 
mass by numerically solving the BA equations (0) and 
0. Our results are displayed in Fig^ Note that when 
solving the BA equations, one should take care of choos- 
ing N~'^ < 7 < TV^. Indeed, if 7 < TV"^, the potential 
energy per particle in the weak coupling limit is lower 
than the zero point kinetic energy fi^/mL^ and the gas 
is therefore non-interacting (effectively 7 = 0). In the 
strong coupling limit and for the same reason, if 7 > N'^, 
the system behaves as a TG gas (effectively 7 = 00). 

We now turn to harmonically trapped bosons in the 
TG regime and rely on spin hydrodynamics introduced 
for uniform systems l^^. As the ground state is fully 
polarized we assume the equilibrium (longitudinal) spin 
density S{x) — n{x)e3 and study small transverse spin 
density fluctuations 6S{x,t) = SSiBi + (55262, where 
61:62,63 form an orthonormal basis in the spin space. 
For a large N, the equilibrium density profile n{x) in a 
harmonic trapping potential V{x) — muP'x^ j2 is given 
by the Thomas-Fermi expression 



n{x) — nQ\/\ — {x/ RY 



(19) 



Here tiq = n,(0) is the density in the center of the trap 
and R = ^2hN jmuj is the Thomas-Fermi radius. For a 
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strong but finite coupling Eq. H19(l represents the leading 
term, with corrections proportional to inverse powers of 
7o — ■mg/h?nQ. The spin density fluctuations SS obey 
the following linearized Landau-Lifshitz equations |2^ : 
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n{x) 6S2.1 
2 TO*(a;j n{x) 



(20) 



In the TG regime the effective mass entering the equation 
of motion (|20|1 depends on the density profile n{x) as 



m*{x)/m « 3j{x)/2tt'^ ~ 3mg/2n^h'^n{x). 



(21) 



Using the density profile (|19|) and introducing a complex 
function 



n{x)^{x,t) = ^5*1(3:;,*) +iSS2{x,t), 



(22) 



one obtains from Eqs. H20|l : 



6 7oA^ Vl - 



dx(l~X^)dx^.. (23) 



where X — x/R is the dimcnsionless coordinate, and 
we assumed the stationary time dependence ^{X,t) — 
e-^^^^{X). Equation shows that the typical fre- 
quency scale of the isospin excitations is given hy lo/^qN, 
which is smaller than the scale uj of acoustic frequencies 
by a large factor 7oiV. The exact solution to this equation 
was obtained numerically using the shooting method, and 
the spectrum shows only a small difference from the semi- 
classical result 



1^ 



j = 0,1,2,, 



(24) 



where the numerical factor is A = 7r5/48r''(3/4) « 2.83. 
For w ~ 100 Hz, 70 10 and N - 100 as in the ex- 
periment 0, the lowest eigenfrequencies f2j are two or 
three orders of magnitude smaller than acoustic frequen- 
cies and are ~ 0.1 Hz. 

In conclusion, we have found extremely slow (iso)spin 
dynamics in the strong coupling TG regime, originat- 
ing from a very large effective mass of spin waves. In 
an experiment with ultra-cold bosons, this should show 
up as a spectacular isospin-density separation: an ini- 
tial wave packet splits into a fast acoustic wave traveling 
at the Fermi velocity and an extremely slow spin wave 
[2^ . One can even think of "freezing" the spin transport, 
which in experiments with two-component ID Bose gases 
will correspond to freezing relative oscillations of the two 
components. 
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